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Abstract 

>: 

O , Recent works proposed the relaying at the MAC layer in cognitive radio networks whereby the primary packets 



are forwarded by the secondary node maintaining an extra queue devoted to the relaying function. However, relaying 
of primary packets may introduce delays on the secondary packets (called secondary delay) and require additional 
power budget in order to forward the primary packets that is especially crucial when the network is deployed using 
sensors with limited power resources. To this end, an admission control can be employed in order to manage efficiently 
the relaying in cognitive radio sensor networks. In this paper, we first analyse and formulate the secondary delay 
t/3 ■ and the required power budget of the secondary sensor node in relation with the acceptance factor that indicates 

whether the primary packets are allowed to be forwarded or not. Having defined the above, we present the tradeoff 
between the secondary delay and the required power budget when the acceptance factor is adapted. In the sequel, we 
^ \ formulate an optimization problem to minimize the secondary delay over the admission control parameter subject to 

a limit on the required power budget plus the constraints related to the stabilities of the individual queues due to their 
interdependencies observed by the analysis. The solution of this problem is provided using iterative decomposition 
methods i.e. dual and primal decompositions using Lagrange multipliers that simplifies the original complicated 
problem resulting in a final equivalent dual problem that includes the initial Karush Kuhn Tucker conditions. Using 
■ the derived equivalent dual problem, we obtain the optimal acceptance factor while in addition we highlight the 

possibilities for extra delay minimization that is provided by relaxing the initial constraints through changing the 
values of the Lagrange multipliers. Finally, we present the behaviour of secondary delay in case infinite and finite 
' queues are employed and thus the overflow and blocking probabilities are assessed respectively. 

Index Terms 

packet relaying, queues, cognitive radio sensor networks, optimization, Lagrange duality, decomposition methods. 



I. Introduction 

Cognitive radio (CR) has attracted a lot of interest in the last decade from research and industrial communities of 
communications and networking. In general, a CR system allows the spectrum utilization of a licensed frequency 
channel by a secondary system (unlicensed) without significantly affecting the efficiency of spectrum utilization of 
the primary (licensed) user [l]. In order to exploit this functionality, a new type of deployment is needed using 
nodes that acquire the knowledge of channels occupancy. The new demands for CR networks (CRNs) deployment 
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have shown via the standardization activities that a new type of wireless sensor networks (WSNs) will be shaped 
targeting to exploit the benefits of the efficient spectrum utilization This type of CRNs were named cognitive 
radio sensor networks (CRSNs) that require a customized handling in order to achieve the desire results O. 

Several issues are raised from such a new application of WSNs such as the error free spectrum sensing (SpSe), 
new communication protocols based on the SpSe results ||4l and new resource allocation mechanisms as packet 
relaying through the secondary nodes {S). Packet relaying has originally proposed first in |6| in which a secondary 
node retains two queues one for its own secondary packets and one for providing the relay of primary packets. 
Based on this concept, the authors studied the maximum allowable throughput achieved by the cognitive user while 
maintaining the stability of the overall distributed queueing system. However, the delay imposed on the secondary 
packets was not studied. Moreover, the additional power budget required for the relaying process were not taken 
into account that is very important in case of CRSNs with limited power resources. A recent work in Q has 
considered a cognitive relaying framework that manage to forward relaying information; however, although delay 
expressions were given, still power constraints were not considered and an optimization solution is not provided in 
general. A more recent work in IS] presents the tradeoff between the sensing and the energy consumption with 
relays; however, the relaying process is not taken place at the MAC layer but at the physical layer looking into 
corresponding cooperative strategies. In |9 j, authors first describe the delay and the power consumption required by 
the relaying process in CRSNs and in this paper we extend this work defining and solving a specific optimization 
problem highlighting useful insights in this new topic. 

Specifically, we analyse the delay of secondary packets and the power consumption which are both requirements 
for the implementation of the packet relaying at the secondary node. We consider that the secondary node (i.e. 
cognitive) retains its own traffic in the secondary queue while the relayed primary packets are stored in a separate 
devoted primary queue. The secondary node is able to control the admission of primary packets passed through 
the relaying primary queue using an admission control parameter. We assume a scheduling strategy that provides 
priority to relayed primary traffic over the secondary one. Based on this cognitive radio model, we obtain the 
formulation of the secondary delay and the required power budget whereby we highlight the following relaying 
trade-off: on one hand, the admission of more primary packets will reduce the secondary delay and on the other 
hand, it will increase the required power budget. This particular tradeoff can be investigated taking into account 
different objectives. In our case, we opt to formulate a minimization problem for the secondary delay over the 
admission control parameter i.e. acceptance factor subject to the constraint in the power budget and the separate 
stabilities' rules of each individual queue. Using optimization theory and decomposition methods in particular, 
since the original is complicated and depends on several constraints, we prove that indeed an optimal acceptance 
factor value exists under several conditions in which the minimum delay is achieved and yet the power budget is 
sufficiently kept in the specific upper bound. Finally, we discuss the relaying behaviour considering different queue 
models for the secondary queue in terms of buffer size (i.e. infinite or finite) giving more light from practical point 
of view. 

The rest of this paper is organized as follows: Section II provides the system model and the necessary assumptions. 
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The analysis is provided in section III, and in section IV we formulate the optimization problem and in the sequel 
we provide its solution. Section V presents the different behaviour of the relaying queue with infinite and finite 
buffers while the simulation results are presented in section VI. This paper is concluded with section VII. 

II. System Model 

We consider the cognitive radio sensor network depicted in Fig. 1 that consists of a pair of a primary user (PU) 
transmitter and receiver denoted as PU-Tx and PU-Rx respectively and a pair of a secondary transmitter (SU-Tx) 
and receiver (SU-Rx). There are three hnks of interest namely the primary Unk between the PU-Tx and PU-Rx 
with a channel gain denoted by gp, the hnk between the SU-Tx and PU-Tx with a channel gain denoted by gs,p 
and the Unk between the PU-Rx and SU-Tx with a channel gain denoted by Qp^s- All channel gains are assumed 
to be zero mean circularly symmetric complex Gaussian random variable with variance and independent for 
all i where i takes the values 'p', 's', 'ps' and 'sp' indicating the primary, secondary, primary to secondary and 
secondary to primary links respectively. Each link is affected by complex additive white Gaussian noise AWGN 
with zero mean and unit variance and independent for all links. We assume that successful transmission is achieved 
when the instantaneous signal to noise ratio (SNR) on the i — th link, employing the maximum transmission power 
Pi, given by | gi p Pi, exceeds a predefined threshold 'yth,i- 

The considered scenario throughout this paper is that the SU-Tx utiUzes the spectrum resource in order either 
to serve its secondary receiver (SU-Rx) whenever the primary hnk is idle or to serve the PU-Rx by relaying the 
primary packets in case the primary link is not able to achieve adequate communication (i.e. outage condition). In 
this fashion, the SU-Tx plays two roles in the presented CR scenario such as giving on one hand the opportunity 
to forward the primary packets whenever the primary hnk experiences outage conditions and on the other hand to 
exploit the spectrum for itself whenever is not used by the primary hnk thus forwarding the secondary packets to 
its own user. Henceforth, we consider that the SU-Tx maintains two queues denoted as Qs and Qps for serving 
its own packer i.e. secondary packets and the primary packets respectively. Moreover, we assume that the PU-Tx 
retains a queue Qp to store its own primary packets. Afterwards, the presented model assumes that the SU-Tx has 
the capabihty to relay or not to relay the primary packets and this is a matter of making a decision taking into 
account the current channel conditions. 

For simplicity we assume that all packets are transmitted in a time-slotted fashion where they have the same 
size and one time slot is required for transmission of a single packet. Initially, we assume that all queues are of 
infinite lengths. The packet arrival processes of self-traffic at each node are independent and with mean arrival rates 
Xp and As (packet/slot) for the Qp and Qs queues respectively. The MAC layer is assumed to obey the following 
protocol. At a certain time slot and given the priority to grant the PU unconditional access to the channel, the PU-Tx 
transmits a packet to the primary destination i.e. PU-Rx. A correctiy received packet by the primary destination 
is acknowledged by sending an ACK message resulting in the primary queue Qp dropping this packet. In case 
of packet failure (indicated by ACK timeout) and assuming that the SU-Tx can hsten to all PU-Tx transmissions, 
the SU-Tx can accept a fraction of the undelivered primary packets in its relaying queue Qps, if it was able to 
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Fig. 1. System model of cognitive radio sensor network with packet relaying capabilities 

correctly decode them and send ACK messages to the primary destination PU-Rx. In case both the SU-Tx and the 
PU-Rx fail to decode the primary data, a retransmission of the packet is initiated by the PU-Tx. We assume that 
the overhead for transmitting the ACK and NACK messages is very small compared to packet sizes. In addition, 
we assume that the ACK messages are always decoded perfectly at the PU-Rx and the SU-Tx. At the beginning 
of each time slot, the SU-Tx is allowed to sense the channel and if it is declared idle, the SU-Tx transmits from 
one of its queues under the conditions discussed below. 

III. Performance Analysis and Metrics of Interest 

In this section, we present, the analysis and the corresponding metrics of interest for the considered cognitive 
sensor node with packet relaying capabilities. We tackle with two performance metrics and obtain their formulation 
and how these metrics are involved in the decision making mechanism at the SU-Tx for deciding on relaying the 
primary packets. The metrics considered are the SU-Tx's delay introduced to the secondary packets in the Qs, 
called secondary delay and the power budget required by the Qps to forward the primary packets. 
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A. Relaying Factor and Scheduling 

We assume a relaying factor / e [0, 1] that is managed by the SU-Tx for making the decision to relay or not 
relay the secondary packets. Increasing / to 1 means that no admission control is accomplished at the Qps thus 
relaying with any restriction is taken place and making / equal to means that no relaying is accomplished thus 
the packets are blocked from the assumed admission control mechanism. Relied on this factor, the SU-Tx performs 
the following scheduling strategy called Prioritized Cognitive Relaying (PCR) where the primary packets have 
always priority against the secondary ones IfTOl . In particular, using this strategy, when an idle slot is detected, 
the SU-Tx attempts to access the primary channel in order to transmit first the primary packets using the relaying 
queue Qps- If possible, the SU-Tx uses the minimum transmission power level sufficient to achieve the target 
(Signal-to-Noise-Ratio) SNR threshold ^th,sp at the primary destination which is equal to ^th,spl \ Qsp Thus, an 
outage probability is defined for each link denoted as Pout,i with i the index for each link. If the link between the 
relay and the primary destination experiences a deep fade or the relaying queue is empty, the SU-Tx will transmit 
from the secondary traffic queue Qs (its own traffic) with the same procedure. According to the adopted PCR, the 
arrival and service rates of the CRSN denoted as Ai and for each i queue respectively encompass the relaying 
factor / in their formulation as follows: 

- Primary Queue Service Rate Hp-. The primary queue Qp will forward the packets once the primary link does 
not experience outage and when the primary to secondary link does not also experience outage taking into 
account the relaying factor and thus it is obtained as 

l^p (1 Pout.p^ ^ f Pout.pi.^ Pout,ps) (1) 

- Relaying Queue Arrival Rate Xps'. The relaying queue Qps will arrive the packets once the primary link does 
not experience outage taking into account the relaying factor and when the primary to secondary link does not 
experience outage taking into account the outage of the primary link and thus it is obtained as 

^ps f^out.p f ^out.pi^ Pout,ps) (2) 

^J,p 

- Relaying Queue Service Rate ^ps'. The relaying queue Qps will forward the packets once the secondary to 
primary link does not experience outage and thus it is obtained as 

IJips = {I - - Pout,ps) (3) 

- Secondary Queue Service Rate ^s'- The secondary queue Qs will forward the secondary packets once the 
secondary to primary link does not experience outage and when the secondary link is also free of outage and 
thus it is obtained as 

= (1 - ^)(1 - PouUs){l - ^(1 - PouUsp)) (4) 



where the outage probability Pout,i for a given threshold in SNR ^th.sp is obtained as follows 

Pout,. = Pr[\ gsp ? P^ < gth,^] = 1 - (5) 
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B. Queuing Modeling and Secondary Delay 

In this section, we analyse the secondary delay i.e. the delay experienced at the queue Qs based on the considered 
traffic models for all queues in the system. Firstly, we will deal with the primary queue Qp. Per the discussion 
given in ifTTIl and lfT2]| about the dominant systems, it is concluded that the primary queue can be uncoupled from 
the overall system of queues. This isolation allows us to model the primary queue as an M/M/1 queue. Hence, 
the expected value of the delay imposed on the primary packets in the primary queue is given by; 

1 - A„ 

Notably, each node chooses its arrival rate based on the application such that satisfying the stability conditions 
of its queue. For more detailed discussion about the stability and its conditions, the reader is referred to lfT2l . 
For the relaying queue Qps and the secondary queue Qs, we model their traffic model by adopting an M/G/1 
priority queuing model. While their inter-arrival times are Poisson random variables, we can model their service 
times iSps and Ss) as Geometric random variables based on the description of the scheduling strategy in |10| . 
The Geometric random variables parameters Pps and ps are given by jips and /is for the relaying queue and the 
secondary queue respectively. The first and the second moments of the average service time for the relaying queue 
{E[Sps\ and E[Spg]) are given as follows 

E[Sps\ - — (7) 

Pps 

E[Sl] = (8) 

Pps 

With the same logic, we can obtain the first and the second moments of the average service time for the secondary 
queue i.e. E[Ss\ and E[Sl]. 

We are now ready to analyse and derive the secondary delay based on the delay experienced at secondary packets 
i.e. secondary queue Qg that can be produced by three distinguished factors according to the model presented in 
ifTSl . Taking also into account the scheduling strategy presented above, the different delay factors for a secondary 
packet of interest inside the secondary queue is assembled as follows: 

1) The waiting time for the currently transmitted packet being served in the secondary queue Qs i.e. once the 
packet of interest is arrived well known as the residual time denoted as Di 

2) The waiting time for being served the packets already appeared in the secondary queue Qs upon the arrival 
of the packet of interest denoted as D2 

3) The waiting time for being served the primary packets either from the primary queue Qp or the relaying 
queue Qps that will be appeared in that queues after the arrival of the packet of interest in the secondary 
queue Qs denoted as Da 

We donate as pi the utilization factor for each queue i that is given in general by the formula p = \iE[Si\. The 
expected value of the delay Di will be the sum of the residual time at the secondary queue and the primary queue 
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obtained as follows 

' 1-Pp 2E[Ss] 1-Pp 2E[Sps] 
where pj/(l — pp) is the conditional probability, that the packet of i — th user is being transmitted given that no 
primary packet is transmitted. We consider now the second component delay D2 which as explained is due to 
other packets in the secondary queue found by the new arrival. Considering the Littles law of queuing theory on 
which the long-term average number of customers in a stable system is equal to the long-term average effective 
arrival rate, A, multiphed by the average time a customer spends in the system, we obtain the following formulation 
considering the summary of the waiting times introduced from the primary and secondary queues respectively 

D2 = PpDp + p,D, (10) 

In the same concept i.e. the Littles law, we can define the third component delay where only the delay from 
the primary packets is considered as follows 

D3 = PpDp (11) 
By adding the three components we can obtain the delay of the secondary queue Qs that eventually become 

1 - Ps 

C. Relaying Power Budget 

Except the delay at the secondary queue Dg, the other important metric for our model is the power consumed 
on the relayed traffic relative to the total power consumed by the SU-Tx. An expression for the average power 
consumed by the SU-Tx on its own traffic as well as the relayed packets is analysed below. To this end, we define 
the relaying power budget as the ratio between the average power per slot consumed on the relaying effort Preiay 
and the total average power per slot consumed Ptotai that is expressed as follows 

p _ Preiay __ Apg£'[Pgp] 

Ptotai \psE\P,p\ + \,E\P,\ 

where E\Pg^ is the average transmission power per relayed packet that is obtained as below considering the 
assumption of the considered CRSN. First, the transmission power consumed on a relayed primary packet is 

where it is a random variable since it depends on the channel gain | ggp P which is a chi-square distributed random 
variable. To calculate the average power per packet consumed on relaying, we need the probability density function 
(PDF) Pgp that is defined as follows 

= ^e^, Va; > (15) 
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where x ~\ gsp P that is obviously an exponential random variable. Considering both (fT3] l and (fl4] l. we can assume 
the following for the PDF 



= VP,p>0 (16) 



, P2 

sp sp 



Since the Psp cannot exceed the maximum limit P*p then the average power can be obtained as follows 

sp 



E[Psp < p,;] = p{x) = / Pspp{p,p)dp,, 

Jo 



"^^e ^^^dP.sp,yPsp>0 (17) 

^ spJ^sp 

Let z = dPsp 2^^^ that will give the lower hmit of the integration as equal to 

upper limit as infinite where finally we have 

J ^ 2^ / z-ig-^dz (18) 



The integral in (fTSl l is known as the upper incomplete gamma function and is calculated using The Exponential 
Integral lfT3l denoted as Ei [.] and thus we have the following formulation 

T?\r> 1 lth,sp I 1th, sp ^ 1th, sp ,,„s 
'^sp '-'sp^sp sp 

Substituting fT% into (T3[ and maintaining some simplification, yields 



p _ ^pslth,sp£sp(^ s ^2Q) 

^ sp 



\slth,spespOl + \slth,sSsOl 



IV. Relaying Tradeoff at MAC Layer of Secondary Node 

In the previous section, we analysed and established the formulation of the secondary delay Z)^ and the required 
power budget F of the secondary queue Qs in relation to the acceptance factor / of the admission control policy 
applied at Q^. In this section, we present the relaying tradeoff by which the decision making to relay or not to 
relay will be built. This decision is based on the fact that on one hand increasing the acceptance factor / the 
average delay increases and on the other hand the power budget increased is proportional to acceptance factor 
as depicted in Fig. 2 for some specific values of the considered network. 

This behaviour led us to the definition of a minimization problem over the acceptance factor / considering the 
average delay as an objective function subject to a constraint on the required power budget F. Taking into 
account the constraints imposed by keeping the stability of each queue, then the optimization problem can be 
defined as follows 



9 




Fig. 2. To relay or not to relay tradeoif based on the average delay Ds and power budget F 



mm Dsif) (21) 

Xp < fip 

S.t. < 

As < ^J-s 
r <Tth 

This problem is highly complicated and thus difficult to solve in this form due to the following reasons: 

1) The objective function Dg is not easily differentiable over the acceptance factor / 

2) There are four complicated coupUng constraints related to the acceptance factor / 

3) The expectation is taking over three random variables gss,gsp and^ps 

Henceforth, we will adopt the concept of decomposition methods for solving this optimization problem ifTSll . We will 
use decomposition methods in order to decompose the original complicated problem into equivalent sub-problems 
as presented in |[T1 IfTSl . The existing decomposition techniques can be classified into primal decomposition and 
dual decomposition methods ifTSl . The primal decomposition provides the decomposition of the original primal 
problem and it is appropriate when the objective function is formulated over two parameters. On the contrary, 
the dual decomposition is based on decomposing the Lagrangian dual problem derived from Lagrange multiplier 
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application ifTSjj . The problem is convex as it is proved in the Appendix although using decomposition the dual 
problem can always be considered convex. 

In our problem, we will follow an hierarchical decomposition that applies to our layered architecture by changing 
the primal and dual decompositions recursively. In particular, the basic decompositions are repeatedly applied to 
the problem to obtain smaller and smaller sub-problems. This approach corresponds to first applying a full dual 
decomposition, and then a primal one on the dual problem. As illustrated in Fig. 3 the master problem is decomposed 
into two levels using dual decomposition for the two pairs of constraints where the first level includes one lower 
and one higher master problems and the second level includes the dual and primal decomposition of the resulted 
master problems. 

A. First level decomposition 

For the problem in ( l2Tl i. we will use first dual decomposition using the following two pairs of constraints. Which 
pair of constraints will be used is a matter of discussion and depends on how the system should be converged and 
stabilized. In general, using decomposition methods, the convergence and stability are guaranteed if the lower level 
problem is solved on a faster time scale than the higher level problem. Thus, we choose to have at the low level 
the stabilities of primary and secondary queues Qp and Qs since we want to converge first in these constraints 
and at the high level the stability of secondary-primary queue Qps and power budget F is retained that are more 
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significantly demanding for the considered optimization problem. 

Going into details, while working at the first level, we obtain the Lagrangian function of ( |2TI ) considering the 
constraints of the stability of secondary-primary queue Qps and power budget F that can be written as follows 

Vi,V2) = i5s + Vi{\ps - fJ-ps) + 1^2(X - Tth) (22) 

in which the Lagrangian dual problem is based and obtained as follows 

mm Ds + i^i^Xps - fips) + i^2i^ -^th) (23) 

{Xp < /ip 
Xs < Us 

The dual function in ( l22b serves as a lower bound on the optimal value of the primal problem ( |20l ). For instance, if 
we denote the minimization problem in ( |23] ) as gi{i'i, 1^2) and the optimal value of ( 1211 1 by c?*, then the inequality 
dl > 1^2) holds for any non-negative values of i^i and 1^2- Having defined the lower master problem as in 

( l23T l then the higher master problem using the dual optimization problem of (I2TI) is formulated as 



max gi(z^i,j/2) (24) 



S.t. 




This approach provides in fact the solution of the dual problem instead of the original primal one and it will only 
give appropriate results if strong duality holds. The problem in (l24l) is convex even if the original problem is not 
convex because it is the point wise maximum of a family of affine function {vi, 1^2), and thus the KarushKuhnTucker 
(KKT) conditions are satisfied and the duality gap is indeed zero |fT6l . Afterwards, we can solve the primal problem 
by solving the dual problem based on the dual function 51(1^1, 1^2)- This type of problem is solved by using iterative 
gradient or sub-gradient algorithms according to the differentiable capabilities of the dual problem. In our case, 
since the dual function is not differentiable, we will choose the sub-gradient solution |18|. The iterative algorithm 
will be described at the end of this section in conjunction with the solution of second level decomposition described 
in the next section. 



B. Second level decomposition 

After the dual decomposition performed at the first level, the derived lower and higher master problems ex- 
pressed by (l23T l and (|24] | respectively will be solved using dual and primal decomposition methods respectively. 
More specifically, the lower master problem presented by (|23T l will be solved using dual decomposition since the 
optimization parameter is still the acceptance factor /. Thus, we obtain the Lagrangian function of ( |23] l considering 
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the constraints of the stability of secondary and primary queues that can be written as follows 

L2 if, ^ Ds+ Vl {Xps - f^ps ) 

+l^2{r-Tth)+aK~^^s) (25) 

in which the Lagrangian dual problem is based and obtained as follows 

mm Ds + vi{Xps - tJ-ps) + i^2{^ - ^th) + ^(As - Ms) (26) 

s.t. < lip 

For the known reasons discussed previously for the first level decomposition, the dual optimization problem of ( l23b 
is written as 

max .92 (C) (27) 

where 52(C) is the minimization problem in ( |26] |. Again, as we will see below the dual problem in (l27T i will be 
solved by using an iterative sub-gradient algorithm. 

The higher master problem at the first level expressed by (l24l) will be solved by using primal decomposition since 
the optimization is over two parameters i.e. vi and V2. This will be separated again into two levels of optimization 
for each optimization parameter and thus two sub-problems are produced named lower and higher sub-problems. 
In the same notion, which sub-problem will be considered lower or higher depends on the convergence conditions 
that we are willing to retain. To this end, we prefer to keep at the higher level the optimization over the Lagrange 
multiplier V2 that is related to the power budget constraint. Afterwards, the lower sub-problem is defined as follows 
over the Lagrange multiplier vi 

max gi{vi,V2) (28) 
s.t. \vi >0 



while in parallel at the higher level, we have the following problem in charge of updating the coupling variable V2 
by solving 

max gl{vi,V2) (29) 

S.t. IV2 >0 



where gi{vi, V2) is the optimal objective value of problem ( l28l l for a given V2. Both lower and higher sub-problems 
in ( |28] ) and ( |29l ) are convex since the original problem in ( |24] | is a convex optimization as represents the dual 
problem of the master primal problem. Therefore, based again on sub-gradient method, we will use an iterative 
algorithm where the lower problem in ( l28T l will converge first for given V2. Details will be given below for the 
overall algorithm. 
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If we gather all aforementioned decompositions in one formula, we will realize that the equivalent optimization 
problem of the original problem in ( |2TI ) can be written as follows 

max max max min Ds + i^iiXps — ^J■ps) 

i'2>0 i^i>0 5>o / 

+Mr-Tth) + a>^s-^^s) (30) 

s.t. < fj.p 

Equation (l30l l substitutes the original optimization problem in (|2T]) where the constraints have been introduced 
multiplied with the Lagrange multipliers 1^1,1^2 and that are actually the KarushKuhnTucker (KKT) conditions and 
they can be considered as shadow prices giving a freedom to potentially improve our objective function i.e. minimize 
the average delay. The additional decrease in the objective function is accomplished due to the relaxation of some 
of the given constraints [,19|. Usually, how much performance improvement we can achieve by over-optimizing 
the objective function is determined by an allowable gap. However, as we will see from the simulation results, we 
can achieve better relaying tradeoffs without breaking the initial rules i.e. the constraints. All this discussion will 
throw more light on the formulated relaying tradeoff by looking into the sensitivity (shadow prices) analysis based 
on the formulated KKT conditions that retains the restrictions of the initial constraints ll20l . Finally, the equivalent 
problem in (|30] | will be solved by the following iterative subgradient-based algorithm as has been derived from the 
hierarchical decomposition described above. 

Algorithm 1 Hierarchical Decomposition 

1) Initialization: and k — 1 iteration 

2) Repeat the initial values 

a) Initialization: and j = 1 iteration 

b) Repeat the initial values 

i) Calculate the acceptance factor / from ( |26] ) 

ii) Calculate the sub-gradient at ^^ (As — /is) 

iii) Update ^■'+^ by = ^-^ + a^-'(As — /is), where a is the step size 

c) Stop once | — ^-^ |< e , where e is the convergence rule 

d) Calculate the parameter i/^ using (29) based on the previous values ^ and / 

e) Update 1^2^^ by v^^^ = i^^ + '^i'^i)^ where a is the step size 

3) Stop once | 1^2^^ ~ ^2 I — ^ ' where e is the convergence rule 



The algorithm above is based on the subgradients of .gi(i^i, 1^2) and 52(0 given by the following propositions. 
Proposition 1: The subgradient of gi{i>i,U2) is si(/) — Hs — for the j — th iteration and Si{f)y^^u2 is an 
element of 9si(j^i, !^2) Q 

^ ds\{v-i,U2) denotes the set of all subgradients at u\ and U2 that is called the subdifferential. 
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Proof: For any (/ii,/i2) G dom{gi) , since gi(/ii,/i2) is obtained by maximizing /ii, /i2) over / G 

dom({D)s), we have gi{fii,ii2) > -^i(/iyi.i/2 : /^i: A'z) S lfT6l ll2n . Moreover, since /,yi,,y2 achieves the maximum, 
we have 1^2) ~ Li{fy^^^^,^i,^2)- Combining the pieces, we obtain 



51(^1,^2) > il(/l.l,<.2>Ml,M2) 
= Li{fi,^^v2iVl,V2) + Li{f^^^„2, 111, fl2) - Ll{fyi,V2iVl,V2) 

= gi{>^i,v2) + H{f^^,^^Y' {{^ii ~ H2) - {vi - V2)) (31) 

Using a similar methodology, we can prove the existence of the subgradient g2{C) ™d then the problem can be 
solved by algorithm 1, which requires the calculation of the subgradients gi{vi, V2) and 52(C) each of their own 
iteration. 

V. Relaying Behaviour for different M /G/l Queuing Models 

One important issue in networks in general is to model the traffics of the primary users and secondary users 
Il22l . In this section, we will translate the acceptance factor / into a specific metric of MjGjX queuing model 
throwing more light on the performance impact in the considered cognitive radio model under realistic applications. 
More specifically, since the acceptance factor determines whether or not a new primary flow can be admitted to the 
secondary queue, we can merge it in the overflow probability and blocking probability for the infinite and finite 
buffer cases respectively [23 1. 

Infinite huffier. Although in case of infinite buffer, the notion of admission control can not be claimed, we can 
adopt the term of overflow or tail probability that is defined as the probability in which an arriving packet finds no 
room in the hypothetical queue. This aspect is identical to the example where in a waiting room, there is limited 
number of chairs for the customers, and any customer that upon arrival does not find a seat has to stand |24|. Lets 
denote the limit of packets equal to K and thus the overflow probability is expressed as follows 

00 K 

Vov = Pr\N>\\= p„ = l-Y,Pn (32) 

n=K+l n=0 

where N is the total number of packets in the system and p„ is the queue state, namely the probability that there 
are packets in the queue. The probability p„ is obtained as the Laplace inverse transform of its generating function 
defined as follows 

= "-f'^-^-'"-' ,33, 

B{\ — \z) — z 

which represents the PoUaczek-Khinchin formula ll25l . where the Laplace inverse transform when Poisson arrival 
process and generalized service time are assumed results in the following expression for the probability pn 



where fv-^.v^ is the value of the acceptance factor / for the Lagrange multiplier values v-i and V2- 
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Substituting ( [34l ) into ( [32] l. we can derive the overflow probability for the infinite M/G/1 queue that is associated 
with the acceptance factor through the arrival rate and the utilization factor of the queues in general. 

Finite buffer. In case of finite buffer the blocking probabiUty plays the role of admission control metric in 
admission control policy. A well-known example of finite systems is the blocking probability in cellular networks 
where the limit resources can be fully allocated. There are numerous approximations for the blocking probability 
possible for M/G/l/K queuing systems as described in 1261 . A good approximation of the blocking probability 

is the one that is based on the overflow probability of an infinite buffer queue that we discuss above and thus 
we can associate it with the analysis of the overflow probability. To this end, the blocking probability is defined as 
follows 

(1 - p)po^ 

Pb = — (35) 

1 - PPov 

where Pov is the overflow probability obtained by ( [32] ) using ( [34] l and eventually the blocking probability is associated 
with the acceptance factor as mentioned above. 

In order to solve the problem in ( |2TI ) over the overflow probability pov defined in ( |32] | and the blocking probabiUty 
Pb defined in [351 we need to discuss for the convexity of these equations over the utilization factor p that is a linear 
function of acceptance factor / and thus we can infer that are convex over / as well. 

Proposition 2: For the range of value p, < p < 1 the first derivative of Pov{p) is increasing and convex on p. 

Proof: Taking the first derivative of overflow probability pov in (13) with respect to p gives: 



Pov(P) = (p_ 1)2 + ^ (36) 

For < p < 1, it is evident that poV (p) > and thus Pov{p) is increasing in p i.e. is convex in p. 
Regarding the convexity of blocking probability in ([32] l, it can be considered as (1 — p)a;/(l — px) that is convex 
increasing to the interval p G [0, 1] with / (0) = (1 — p) and /(I) = 1 as discussed in ll27l . 

VI. Simulation Results 
We provide below the simulation results of the optimization problem above with the corresponding discussion. 
We assume the following values for the SNR thresholds ^th.p = lth,ps — 1th, sp — 1th, s — OdB and the following 
channel variables <jp ~ AdB^cr^^ — 12dB,<7^p — 8dB and cr^ = 12dB. The transmission power needed for each 
queue is Pp = 1, = 1 and Psp = 0.25 . Fig. 4a depicts the optimal acceptance factor /* and the corresponding 
average delay Ds{f*) vs. arrival rate at the primary Queue Xp for Tth = 0.2 power budget threshold and As = 0.1 
as secondary arrival rate. The figure gives insights of how much relaying is needed for different values of the 
primary loading. For small values of Xp, there is low traffic in the primary queue, and the SU-Tx (i.e. cognitive 
node) has more opportunities to access the channel, so it can accept all the undelivered primary packets for relaying 
while wasting a small percentage of its power As increasing the primary load increases the number of relayed 
packets, the cognitive node reduces the accepted fraction of these packets to keep the predefined power budget. For 
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Fig. 4. Optimal acceptance factor /* and corresponding average delay Ds (/* ) vs. arrival rate at the primary Queue Ap 



high values of Ap, the SU-Tx might choose not to relay the primary packets, as a small reduction in the secondary 
delay is associate with very high power consumption. In the same notion. Fig. 4b depicts the optimal acceptance 
factor /* and the corresponding average delay Ds{f*) vs. arrival rate at the secondary Queue As for Tth — 0.2 
power budget threshold and Ap = 0.3 as secondary arrival rate. The insight is opposite of the previous one i.e. vs. 
arrival rate at the primary queue Ap. For high values of A^ the cognitive i.e. secondary node does not has enough 
opportunities to access the channel since the primary packet are mainly forwarded as pointed out above. 

Fig. 4c depicts the optimal acceptance factor /* and corresponding average delay Ds{f*) vs. power budget 
threshold Tth- The effect of the relaying power budget threshold Tth on the secondary delay Dg is shown. The 
solution for the optimization problem for different values of Tth is plotted. The simulation parameters are the 
same as that for the previous figures and for arrival rates we have obtained Ap — 0.5 and A^ ~ 0.1. The general 
interpretation of this results is that if low power budget is assigned to the relaying channel (i.e. < Tth < 1) , the 
cognitive sensor node accepts a little fraction of the undelivered primary packets (/ < 0.1) which does not result 
in a significant reduction in the secondary delay. However, increasing the budget threshold allows more packets to 
be accepted for relaying and, in turn, reduces the secondary delay but with a limit. The figure shows that, spending 
more than 0.45 of the total power on relaying, will not result in any reduction in the secondary node delay as it 
accepts all the undelivered primary packets for relaying. 

Fig. 5a depicts the sensitivity of the average secondary delay Ds in terms of the KKT conditions of the original 
optimization problem in (l2Tl i as they introduced by the dual problem in (|30] |. We keep the same values for the 
parameters as previously. In three axes the KKT conditions i^i(Aps — fips), '^2{T — Tth) and ^(A^ — /is) are denoted 
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Fig. 5. Optimal acceptance factor /* and corresponding average delay Dsif) vs. arrival rate at the secondary Queue As 




Fig. 6. Optimal acceptance factor /* and corresponding average delay Z?s(/*) vs. power budget thresliold ^^^^ 
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Fig. 7. Arrivals rate Xp sensitivity of primary queue Qp on shadow prices i.e. Lagrange multipliers 

as vi{f), 1^2 (/) and ^(/) respectively at x, y and z axes as well. Hence, they can be called shadow prices since they 
can provide decision-making with powerful insights into our problem. For instance, it is obvious from the Fig. 5a 
that the KKT condition related to the power budget threshold i.e. 1^2 (/) might not be kept as long as the Lagrange 
multiplier V2 is increased more than a threshold since obviously the V2{f) is increased in this case beyond the zero 
value. This is applied when the primary arrival rate Ap is getting higher i.e Ap > 0.3. In the contrary, the more 
increase of the Lagrange multipliers V2 and ^, the more the decrease of the KKT conditions vi{f ) and ^(/) is 
become and as a consequence the average delay Ds of the dual problem in (|30] l. In Fig. 5b, we depict the KKT 
conditions changing the secondary arrival rate A^ keeping the primary arrival rate Ap ~ 0.1 and we highlight that 
all shadow prices vi{f),V2{f) and f (/) can potentially decrease the average delay without having a harmful impact 
on the stability of secondary queue. 

Fig. 5c depicts the sensitivity of the powers budget F sensitivity on shadow prices i.e. Lagrange multipliers for 
different threshold values Tth- We retain the aforementioned parameters for the CRSN and for the arrival rates 
we assume As = 0.1 and Ap — 0.5. Notably, as long as the power budget threshold Vth is low, the increase in 
i^2(r — Tt/i) KKT condition will not pass the intolerable limit of increasing the average delay. However, the decrease 
is retained for low values of the corresponding shadow price V2{f) and not for the high power budget threshold 
values i.e. Vth — 0.4 and Vth — 0.6. 

Fig. 6 depicts in the top axis the overflow and blocking probabiUties denoted as pov and pb respectively for 
different values of primary arrival rate Ap. We assume that the secondary arrival rate is As — 0.1 and the threshold 
on the powers budget Vth = 0.2. In the same concept, the middle axis depicts the overflow and blocking probabilities 
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Fig. 8. Arrivals rate As sensitivity of secondary queue Qs on shadow prices i.e. Lagrange multipliers 




Fig. 9. Powers budget F sensitivity on shadow prices i.e. Lagrange multipliers 
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Fig. 10. Overflow and blocking probabilities for the optimal acceptance factor vs. primary arrival rate \p, secondary arrival rate As and power 
budget threshold Ttf^ respectively 

Pov and respectively for different values of secondary arrival rate A^. We assume that the primary arrival rate 
is Ap — 0.3 and the threshold on the powers budget Tth = 0.2. Finally, the bottom axis depicts the overflow 
and blocking probabilities pov and pi, respectively for different values of threshold on the powers budget Tth with 
Ap = 0.5 and As = 0.1. From these three figures is inferred that the overflow probability is always more than the 
blocking probability in general and this is due to the relaxing behaviour of infinite buffers against the blocking 
probabihty of finite buffers. On the other hand, the convex properties of both probabilities are confirmed as it is 
highlighted in (26). More remarks on the behaviour of relaying queue are that first the overflow and blocking 
probabihties are not increased up to one level over secondary arrival rate due to the forwarding mainly of the 
primary packets process as pointed out above and finally the probabilities tends to one as long as the power budget 
threshold is relaxed. 

VII. Summary 

In this paper, we consider a sensor system in which a node is able to relay the primary packets at MAC layer 
when the primary cannot achieve adequate communication. The relaying model is based on queue modelhng where 
the stability rules should be satisfied accompanying with the customized rules of our investigation, which are 
the minimum delay induced in the secondary packets and the power budget required for relaying the primary 
packets. The relaying capability is controlled by an admission control factor that is associated with all available 
distributed queues in the system. Based on the considered model, we define an optimization problem with four 
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several constraints. Specifically, the problem is formulated as the minimization of average delay at secondary packets 
as long as the queues stabilities and the power limit are guaranteed. We solve this complex optimization problem 
using hierarchical decomposition method thereby the several constraints are imposed separately applied first and 
second level decomposition, which finally are broken into more concrete sub-problems. This method results in 
the equivalent problem that encompass all the constraint requirements as KKT conditions. We obtain results that 
highhght useful outcomes in terms of system performance as well as indicate the sensitivity strength of the objective 
function both to the separate numerous constraints. We also assess our solution considering infinite and finite buffer 
capacity for the queue that forwards the primary packets. 

Appendix 

Proposition 3: For the range of value /, < / < 1 the first derivative of Dg (/) is negative and increasing (i.e. 
Dsif) < 0) and thus the Ds{f) is decreasing and convex on /. 

Proof: Looking into the average secondary delay Ds{f) in (13), both pp and Ps utilization factors are constants 
while the pps and the average primary delay Dp do not. Hence, differentiating Ds{f) with respect to / gives: 



D,{f) = AppM) + BDp if) + C (37) 

where A, B and C are some constants while Ppsif) and Dp (/) are the first derivatives of utihzation factor and 
primary delay respectively which are obtained as follows 

(1 ^p)Pout,p(,^ Pout,ps) 
(1 Pout^p ~t~ f Pout^piX Pout^ps^y^ 



7=; / e\ \^ ^p )J- out,p\i- J-out,ps) ,r,a\ 

yj) = n — — , fu n ^^^^ 



and 



(1 Pout,p) f Pout,p{^ Pout^ps) ^p 



f^pPout,p{^ Pout,ps)^ ^2^^ 



(1 Pout,p ~\~ f Pout,p{^ Pout,ps}}'^ 

For < / < 1 it is clear that D^ (/) < i.e. negative and increasing as well and thus Ds{f) is decreasing as 
depicted in Fig.2 and convex as a conclusion. 
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